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Ginzburg-Landau Equations for Coexistent States of Superconductivity 
and Antiferromagnetism in t - J model 
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Department of Physics, Kobe University, Kobe 657-8501, Japan 

Ginzburg-Landau (GL) equations for the coexistent state of superconductivity and antiferro- 
magnetism are derived microscopically from the t — J model with extended transfer integrals. 
GL equations and the GL free energy, which are obtained based on the slave-boson mean- 
field approximation, reflect the electronic structure of the microscopic model, especially the 
evolution of the Fermi surface due to the change of the doping rate. Thus they are suitable for 
studying the material dependence of the coexistent states in high-T c cuprate superconductors. 
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1. Introduction 

The discovery of the coexistence of antiferromagnetism (AF) and superconductivity (SC) 
in multilayer high-r f cuprates has stimulated wide interest. 12) Antiferromagnetic superex- 
change interactions in high-T c cuprate superconductors, which are strongly correlated elec- 
tron systems, are thought to be the origin of two ordered states; thus understanding the con- 
dition for coexistence may give insight into the mechanism of superconductivity. 

In single-layer and bilayer cuprates such as La- and Y-based compounds, it has been well 
known that AF is easily suppressed by a tiny amount of carrier doping. 3 4) On the contrary 
in multilayer systems (in this paper the term "multilayer" will refer to three or more layers 
in a unit cell) such as HgBa2Ca4Cu50i2 +v , AF survives up to much higher doping rate and 
coexists with SC state. NMR measurements revealed that the coexistence was not due to a 
proximity effect but a genuine phase transition within a Cu02 plane. 12) Multilayer cuprates 
have flat Cu0 2 planes with a perfect square lattice and are known to be free from disorder in 
contrast to La- and Y-based compounds. Combined with their high T c of more than 100K, 5) 
multilayer cuprates can be viewed as ideal systems to study the mechanism of high T c . In this 
sense it is desirable to explore the nature of the coexistent state of AF and SC theoretically. 
Low-energy electronic states of high-r c cuprates are described by the t - J model. 6_8) 
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In the case of single-layer and bilayer systems the AF order is easily destabilized by strong 
fluctuations due to low dimensionality. Assuming the absence of AF order, mean-field (MF) 
theories 9 ' 10) based on the slave-boson (SB) scheme 1 1) to treat the condition of no double 
occupancy and the gauge theory, 8 ' 12) which takes into account the low-energy fluctuations 
around mean-fields, capture many important properties of single-layer and bilayer high-T c 
cuprates. In multilayer systems, on the other hand, relatively strong three dimensionality 
may stabilize AF order. 13 ' 14) This situation can be suitably treated by MF theories for the 
t — J model by taking AF order into account. Actually MF calculations for the t — J model 
predicted that AF survives up to 8 < 0.1 - 0.15 (6 being the doping rate) and it may coexist 
with SC, 13 ~ 16) and a similar result was obtained by the variational Monte Carlo method. 17 - 1 

In this paper, we derive GL equations and the GL free energy microscopically from the 
two-dimensional t — J model with extended transfer integrals (extended t — J model) based 
on the SBMF approximation. In the MF approach the phase diagram will not be sensitive 
to the number of layers. It is the shape of the Fermi surface, in particular, the condition 
for the nesting that is crucial to determine the occurrence of the coexistent state, and an 
electronlike Fermi surface can lead to the experimentally observed phase diagram. 13 ' 14) In 
multilayer cuprates we expect that such an electronlike Fermi surface may be stabilized as 
one of the Fermi surface due to strong hybridization between layers. This is the reason why 
we treat a single-layer (single-band) model, and we simulate the difference of the Fermi 
surface by including the extended transfer integrals. 

The derived GL theory can be used to investigate the spatial dependence of the AF and SC 
order parameters (OPs) in high-r c cuprates, and it may provide information on the electronic 
states in these systems. For example, near the surface or impurity the OPs are suppressed, 
and their recovery to the bulk values will provide the coherence length, which reflect the 
underlying electronic structures of each system. 

Although the GL theory is reliable only qualitatively except near T c , it can give a sim- 
ple and intuitive description of the coexistence and competition of multiple OPs. Thus, it is 
complementary to more sophisticated methods such as the Bogoliubov-de Gennes and qua- 
siclassical Green's function theory. Previously various models have been employed to derive 
GL equations microscopically; a continuum 18) and tight-binding model 19) with s- and J-wave 
SCOPs, Hubbard model with nearest-neighbor attractive interactions, 20) a model with a spin 
generalized BCS term and Heisenberg exchange term, 21) and the t - J model (without taking 
AF order into account). 22) The method of deriving GL equations in this work is based on that 
by Gor'kov 23) with the extension to include AF order. 2()) 
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This paper is organized as follows. In §2, we present the model and treat it by the SBMF 
approximation. GL equations and the GL free energy are derived in §3. Section 4 is devoted 
to summary and discussion. 

2. Model and Mean-Field Approximation 

We consider the extended t — J model on a square lattice whose Hamiltonian is given as 

H=~Yj tjte^cpte + jYjSj- S f , (1) 

jAo- (j,0 

where the transfer integrals tjt are finite for the first- (t), second- (f), and third-nearest- 
neighbor bonds (t"), and vanish otherwise. /(> 0) is the antiferromagnetic superexchange 
interaction and (j,i) denotes the nearest-neighbor bonds. The magnetic field is taken into 
account using the Peierls phase ip^ = ^ J. A • dl, with A and O = §f being the vector poten- 
tial and flux quantum, respectively. £j a is the electron operator in Fock space without double 
occupancy, and we treat this condition using the SB method 1 !) by writing c ]cr = bjfj a under 
the local constraint Yio-fjo-fja + b\bj = 1 at every j site. Here fj a (bj) is a fermion (boson) 
operator that carries spin cr (charge e); the fermions (bosons) are frequently referred to as 
spinons (holons). The spin operator is expressed as Sj = \ Yia$ fLcapJ 

We decouple Hamiltonian eq. (1) in the following manner. 13 ~ 16) The bond order parame- 
ters (2cr fjg-fccr) and (b]b e ) are introduced and we denote xj<e/2 = (f^fy) = </J/a> for th e 
nearest-neighbor bond. Although the bosons are not condensed in purely two-dimensional 
systems at finite temperature (T), they are almost condensed at low T and for finite carrier 
doping 6(> 0.02). 15) Hence we approximate (bj) « and (b]bj) « 5. The magnetization is 
defined by mj = ^{fj^ff[ - fj±fji), and the superconducting OP on the bond (j, £) (under the 
assumption of the Bose condensation of holons) is given as = {ff\fe,i)- 

In the following we are mainly interested in a region around the tetracritical point where 
the four states, AF, d x 2_ y 2 -wave SC, their coexistence, and the normal states become identical. 
The onset temperature of the bond OPs is much higher than that for AF (T N ) and SC (T c ) in 
this doping region, so that they are almost independent of temperature near the tetracritical 
point. We consider only the spatial variations of rti; and Aj ^ assuming that Xj,i is uniform in 
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space. (Hereafter we denote it as^.) Then the mean-field Hamiltonian is given as 

H MFA = "£[ Z V&W + YXifUrfjr + 'S Z *™fi«J* 
j,cr S=±x,±y S=±x±y 

+ t"8 ''"^^h'la^Tl Z m Mf'nfn - flf'i) 

6=±2x,±2y j,cr j S=±x,±y 

+ ^Z Z i^(fUA-\fUf^ +h - c i +E °> 

j S=±x,±y 

(2) 

with 

E = -Jj] mjm e + J £ (^A/,A;. / + ^| 2 ). (3) 

First we solve the self-consistency equations for x an d me chemical potential p in the 
absence of m, A, and A. Self-consistency equations that determine ;^ and jj. as functions of T 
and 8 are given as 

x = ^ Z (cos pa + cos 5 = 1 - jj Z (4) 

p p 
where ^ p = -{2t8 + ^^)(cos /?. Y + cos /? y )-4f'(5cos p. v cos p y -2t"8(cos 2/^ + cos 2p y )-fi, with 
/ and Af being the Fermi function and the total number of lattice sites, respectively. (Lattice 
constant is taken to be unity.) In the next section we will carry out the GL expansion to obtain 
GL equations for m and A. 

3. GL Equations and GL Free Energy 

In this section we derive GL equations and the GL free energy. The procedure is essen- 
tially the same as that used in ref.20. Coupled equations for Green's functions Go-(j, i, f) = 
-(TrfjAr)/^) and F^O', £, r) = -{T T f j(T (j)f (o .,) can be derived from their equations of mo- 
tion (Gor'kov equations) as 
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J - 

G T (j, t, is„) = G (J, i£n) + 2 Zj Go ^ k ' lSn>> 

x[(A k+Suk + -A kMS ^F\^k + 8i,tJs n ) + m k+6x Gi(k,£Js n )\, 

Gj.0', £ = G (j, i, lE n ) - 2 /j ^ lS n) 

k,Si 

x\(A kM6i + -Afc +5l ^)F^(fc + 5i,€, is n ) + m k+Sl Gi(k, C, is n )\ 
j _ J (5) 



j 1 
F* n (j, £, is n ) = - 2^ G (k, j, -is n ) 

k,Si 

x[{K +Su k + 2 A lk + 6, ) G i( k + 8\,t, ten) + 



2 

where the summation on 8\ (k) is over +x and +y (all sites). Here, Go(j,£,iu)„) is Green's 
function for the system without A and m but with A. GoO', I, ia> n ) is related to Green's func- 
tion for the system without A, Go, as Go(j, I, ie„) ~ G (j, i, is n )e^ 1 , with Go(j, i, ie n ) being 
the Fourier transform of Go(p, is n ) = I /(is,, - ij p ). In the expression of % p , the bond order 
parameter^- and the chemical potential /u determined by eq.(4) are substituted. 

Spin-singlet and spin-triplet SCOPs on the bond (j, j+rf) are expressed in terms of Green's 
functions F^ and F^, 

(Af\j)T = i^jw - fjjj^r = \{Aj, M + A j+nJ y 

(A^u)r = \(f n f j+vl + fjjj^y = \{a^ - A j+Vj y 

= -J^J [ F T^ J + T] ' J ' ien) + + ^ J ' ien) \ 



(6) 



and the staggered magnetization Mj = mp 1 ^'*) (Q = (n, n)) is similarly given using G T and 
Gi, 

Mj = \{f n f n - .r L .i)y 0i: - 

T * (7) 

= 2 ^ G ^' j ' ien) ~ G ^ J ' J ' is "^ e 

We substitute eq. (5) into eqs. (6) and (7) iteratively and keep the terms up to the third order in 
OPs. In the coexistent state of AF and SC, spin-triplet SCOPs that oscillate in a similar man- 
ner as the staggered magnetization may occur, 24 ~ 28) and they are called the 7r-triplet SCOPs. 
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The SCOPs of each symmetry, A s (s-wave), A d (d-wave), and Ar^ (^--triplet /?x(y)-wave), 
can be constructed by making a linear combination of eq.(6), 

A ^) = \ E A?>C/). A d (j) = I[ £ Af (j) - £ Af\j)\ 

ri=±x,±y rj=±x T]=±y (8) 

*So') = 5^0-) + a^o)]. 

Assuming that the SCOPs and M are slowly varying, we take a continuum limit. The 
OPs in the linear terms are expanded in powers of derivatives up tp the second order, and the 
Peierls phase is also expanded in powers of A to the same order. Then after straightforward 
but lengthy calculations we get the following GL equations: 



a s A s + 2&|A/A, - K S (D 2 X + D 2 y )A s - K ds (D 2 x - D 2 y )A d 
+ y x |A^| 2 A, + 2 r2 A>: + y 3 (l A p ,| 2 + |A ro ,| 2 )A. s + 2y 5 (A 2 , + A^A; 
+ y 7 (\A px \ 2 - \A py \ 2 )A d + y & (A 2 px - A 2 py )A* d + y 9 {A px A py + c.c.)A s 
+ 2y n A px A P yA* s + y ms M 2 A s + y spm M(A px + A pv ) = 0, 



a d A d + 2fi d \A d \ 2 A d - K d (D 2 x + D 2 )A d - K ds (D 2 x - D 2 )A S 
+ ri |A,| 2 A rf + 2 72 A 2 A* + y 4 (l A p ,| 2 + \A py \ 2 )A d + 2y 6 (A 2 px + A 2 y )A* d 
+ y 7 (\A px \ 2 - \A py \ 2 )A s + y s (A 2 px - A 2 py )A* s + y w (A* px A py + c.c.)A d 
+ 2y n A px A py A d + y md M 2 A d + y dpm M(A px - A py ) = 0, 



& P \A pX (y) + 0> p 2A py ( x) + 2fl p \A pX (y)\ A pX (y) 

- K p iD 2 x(y) A px(y) - K p2 D 2 {x) A px(y) - K p3 (D 2 x + D 2 )A py{x) 

+ 7 P i \\y(x)\ 2 \x(y) + 2y p2 A 2 py(x) A pxiy) 

+ y P 3(2\A px(y) | 2 A py(x) + A 2 px(y) A py{x) + \A py(x) \ 2 A py{x) ) 



+ (y 3 \A s \ 2 + y 4 \A d \ 2 )A px(y) + 2(y 5 A 2 + y 6 A 2 d )A* 



(9) 



(10) 



(11) 



+ 7v(A,.A* + c.c.)A px{y) + 2y 8 A s A d & PX(y) + (y 9 |A/ + y l0 \A d \ 2 )A py(x) 
+ (TiiA 2 + 7i 2 A 2 )A; vW + (y mpl A px(y) + y mpl A py{x) )M 2 
+ {jspmK ± y dp „Ad)M = 0, 

a m M + 2/3 m M 3 -K m (V 2 x + V 2 )M 
+ (y ms \K\ 2 + y md \A d \ 2 )M + [y mpl (A px \ 2 + \A py \ 2 ) + y mp2 (A px A* py + c.c.)]M (12) 
+ ^JspmVKi A P .v + A^) + c.c] + -y dpm [A* d (A px - A py ) + c.c] = 0, 
where the coefficients appearing in eqs. (9)-(12) are given in the Appendix, and D is the 
gauge-invariant gradient defined as D = V + ^pA. Equations (9)-(12) are the coupled equa- 
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tions that determine SCOPs and the staggered magnetization self-consistently. 

The GL free energy F up to the fourth order in OPs can be obtained from the above GL 
equations in such a way that the variations of F with respect to OPs reproduce eqs. (9)-(12). 
The results are written as follows: 

F = F s + F T + F S t + F M + F SM + F TM + F STM , 
F s = j d 2 r[a s \A s \ 2 + J3 S \A/ + K s 0A s \ 2 + a d \A d \ 2 + A,|A/ + K d 0A d \ 2 

+71 1 A/| A/ + y 2 (A2( At.) 2 + c.c.) 
+K ds ((D x Ad)(D x A s y ~ (D y A d )(D y A s y + c.c.)], 
F T = j d 2 r[a pl (\A% T) \ 2 + |Ajf | 2 ) + a p2 (Ajf (A<f )* + c.c) +^(|A p ./ + |A pv | 4 ) 

+y pl \A% r f\A^\ 2 + y p2 ((A^) 2 (A^) 2 + c.c.) 
+y P 3(\A<*P\ 2 + \A^\ 2 )(A% T \A^T+c.c.) 
+K pl (\D x A% T) \ 2 + \D y A ( pP\ 2 ) + K p2 (\D y AfP\ 2 + |D,Ajf | 2 ) 
+K p3 ((D x A i p n PT(D x A% T) ) + iD y A%Py(p y A<*P) + c.c.)}, 
F ST = J d 2 r[(\A% T f + |A£f | 2 )(y 3 |A,| 2 + y 4 \A d \ 2 ) 

+{((A^) 2 + (A^) 2 )(y 5 (A:) 2 + y 6 (A*) 2 ) + c.c] 

+y 7 (|A$f f - \A ( p n p\ 2 )(A* s A d + c.c.) + y 8 {((A^) 2 - (A<f ) 2 )A*,A* + c.c.}} 
+((A£ r) r Ajf + c.c.)(y 9 \A s \ 2 + y 1() |A/) 
+{Agf A<f (ynCAp 2 + y 12 (A*) 2 ) + c.c.}}, 
F M = J d 2 r[a m M 2 + fr n M 4 + K m (VM) 2 ] , 

F SM = J d 2 r(y ms M 2 \A s \ 2 + y md M 2 \A d \ 2 ), 

F TM = J d 2 r[y mpl M 2 (\AfP\ 2 + |A$ff ) 

+y mp2 M 2 (A% T \A^T + c.c.)}, 
F STM = J d 2 r[y spm MA s (A% T) + Ajf )* 

+y dpm MA d (A^P-A% T ^+c.c.}. 

(13) 

Here, Fs, F T , and F M are the free energy for the singlet and ^--triplet SCOPs, and the stag- 
gered magnetization, respectively, while F ST , F SM , F TM , and F STM describe their couplings. 
Note that F is invariant under all the symmetry operations of the square lattice. F s M and F TM 
are the usual terms to represent the competition of SCOPs and M. F STM is a cubic term that 
couples spin-singlet SCOPs, staggered magnetization, and ^--triplet SCOPS, and it induces 
^■-triplet SCOPs in the coexistent state of AF and SC. Generally in the coexistent state of 
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ferromagnetism and spin-singlet SC state, spin-triplet SCOPs may occur when OPs are not 
uniform in space. 29 33) In the GL theory this can be explained by a cubic term that has a 
gradient coupling of spin-singlet, triplet SCOPs, and the magnetization m. 2()) In the AF state 
magnetization m is oscillating (though the staggered magnetization M is uniform) even in a 
uniform case, and thus ^--triplet SCOP can arise irrespective of the spatial dependence of 
OPs. 

The important point of the present results is that the coefficients appearing in GL equa- 
tions and the GL free energy are determined microscopically. These values depend on the 
parameters of the microscopic model and they reflect the evolution of the shape of the Fermi 
surface. This property can be used to study the material dependence of the coexistent states 
in various multilayer high-r c cuprates. 

4. Summary and Discussion 

We have derived GL equations and the GL free energy microscopically from the extended 
t — J model using the slave-boson mean-field approximation. The derived GL theory can be 
used to investigate the spatial dependence of the AF and SC order parameters in high-r c 
cuprate superconductors. By analyzing the spatial variations of order parameters using the 
present results, information on the electronic states of high-T c cuprates may be extracted. 

A typical example to be studied is the state near the surface or impurity. The interface 
states of heterostructures composed of cuprate superconductors and magnetic materials are 
also worth studying. There the coexistence and competition of superconductivity and mag- 
netism can occur in various ways depending on the materials used. 

Numerical study of the GL equations for the above situations assuming various band 
structure (by choosing the extended transfer integrals) may be interesting, and this problem 
will be examined separately. 
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Appendix: Coefficients in GL Equations and GL Free Energy 

The coefficients appearing in GL equations [eqs.(9)-(12)] and the GL free energy [eq.(13)] 
are given as follows: 

a m = 37(l - — 2j h(p)a) 2 m ), 



817 4 v~i 4 

P 

817 4 v - 1 2? 1 

ri = "8Ar2j /2(p) ^' ^ = 471- 

p 

9 J 2 sr~\ /d€ P \ 2 2 
9J 2 , , ./5f P \2 



p 



SN^ h(P) WJ ***** 

a P i = - J -(l + ^ J^hip) cos 2 p x ), 

p 

J 2 

ap2 = ~AN 2-i 13 COSPx cos Py> 

p 

jA 

P p = 3^vZj /4(/7)cos4 ^' 

p 

jA j 

7pl = W C0$2 Px COs2 Pv ' 7p2 = 4 7pl ' 

p 

jA 

7 P 3 = COs3 Rr COS Py ' ( A>1 ) 

^ = -w^ h{p) (~df> 2qos ' lpx(yh 

p 

J 2 ^ (d%p\ 2 
K » 3 = ~SN Z^ h(p) \dp~J cos ^ cos ^" 



97^ \~\ 2 2 
73(4) = £tt 2^ hip)^^ cos p x , 



p 

A 



SN 



p 



9J 4 2 2 

75(6) = 3^ 2j 7 6^H(rf) COS "^' 



97 4 2 

77 = ^T7 2j h(P)(O s (U d COS 



p 



8N 



p 



974 V r ^ 2 
7s = 2_u h(p)u s (o d cos p x , 

p 

9J 4 v- 2 
79(10) = ^7 l^tsiP^^cos p x cos p y , 

p 

tA 



97 4 2 
711(12) = 2j /6 ^ COSftCOS;,} 

p 

a m = 27(l + — ^/ 7 (p)), 



P lOjTJ 
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6J V-l 

{jms, 7md) = —jT- 2_P<Hp) + Hp)]{0J s , 0J d }, (A-2) 

P 

jA 

{jmpU Jmpl) = ~— 2_j[2I W (p) + 4(P)]{C0S 2 p x , COS p x COS p y }, 

P 

3J 3 \ 

{jspm, Jdpm) = ~— 2_j Hip) C0S P*{ W f 
P 

where a> s = cos /j^ + cos p y and = cos p x - cos p y , and the summation on p is taken over 
the first Brillouin zone. The functions appearing in the integrands are defined as 



Hp) = 


T 2^ G Q (p, is n )G (p, -is n ), 


hip) = 


T 2^ G 2 (p, ie n )Gl(p, -is n ), 


Hp) = 


T 2j G (p, -is n )G (p + Q, is n ), 


Hp) = 


T ^ G 2 (p, is n )Gl(p + Q, -is n ), 


Hp) = 


T ^ Gl(p, is„)G (p, -is n )G (p + Q, -is n ), 


Hp) = 


T 2j G Q (p, is n )G Q (p, -is n )G (p + Q, is n )G (p + Q, -is n ), 


Hp) = 


T 2^ G (p, is n )G (p + Q, is n ), 


Hp) = 


T 2^ G 2 (p, is n )Gl(p + Q, is n ), 


Hp) = 


T 2^ Gl(p, is n )G (p, -is n )G (p + Q, is n ), 


ho(p) = 


T 2j G 2 (p, ie n )G Q (p + Q, is n )G (p + Q, -ie n ), 


I nip) = 


T 2j G (p, ie n )G (p, -is n )G (p + Q, is n ). 



(A-3) 
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